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In this paper a number of explicit a priori inequalities for solutions of the plate equation are 
derived. These inequalities together with the Rayleigh-Ritz technique may be used to compute error 
bounds in various mixed boundary value problems for elastic plates. 



1. Introduction 

In two recent papers [2,3] ' the authors presented methods for obtaining pointwise bounds 
in the three most common boundary value problems for elastic plates. These bounds were of a 
priori type, that is they held for a class of functions required to satisfy only smoothness condi- 
tions. Hence one could approximate the (unknown) solution of one of these problems in terms of 
essentially arbitrary functions, and the inequalities gave bounds on the error. 

In this paper we derive similar a priori bounds in the three most common mixed boundary 
value problems for elastic plates. For simplicity we consider only the case of a simply connected 
region R whose boundary 2 consists of two disjoint portions Si, and S2 (each connected) on which 
different sets of boundary conditions are imposed. It will be clear how the results are to be ex- 
tended if Si, and/or S2 are not connected or if R is multiply connected. 

In this paper we shall restrict our attention to the problem of obtaining bounds for the Li 
integrals of an arbitrary sufficiently smooth function w in terms of L 2 integrals of quantities which 
are data whenever the arbitrary function w is actually the solution a to the problem in question. 
By use of mean value inequalities and the Rayleigh-Ritz technique as indicated in [2,3], the desired 
pointwise bounds are then obtained. The well known Rayleigh-Ritz technique consists in choosing 

N 

w —u— V Oi<£>{, where the <&i are TV linearly independent sufficiently smooth functions, and the a\ 

are determined in such a way as to minimize the terms involving the data of u. 
The particular problems treated here are the following: 

PROBLEM I: A 2 u prescribed in R 

du . ^ 

u, — given on Zi 
an 

u, M{u) given on S2 

PROBLEM II: A 2 w prescribed in R 

du . _-, 

u, — given on Z, 
an 

M(u), V{u) given on £2 



1 Figures in brackets indicate the literature references at the end of this paper. 

*University of Maryland, College Park, Md.: part-time worker at National Bureau of Standards Washington, D.C This research was supported in part by the 
National Science Foundation under Grant NSF-GP-3. 
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Problem III: A 2 a prescribed in R 

u, M{u) given on 2i 

M(u) 9 V(u) given on 22- 

Here A 2 is the biharmonic operator — +2 2 2 + II4 (A denotes the Laplace operator), and 

t— denotes the normal derivative directed outward on 2. Also 
on 

and 

F(„) = -(A u) + (l- ( r)^--(--)j (1.2) 

where — denotes the derivative with respect to arc length on 2, cr is Poisson's ratio and p is the 

radius of curvature. The quantities M(u) and V{u) are proportional to the normal moment and the 
reaction normal to the plate on 2, respectively. 

In subsequent sections we shall refer to the bounding of integrals of an arbitrary function w 
in terms of L 2 integrals of its "data." For example, in problem I the data of w are the quantities 

A 2 w in R, w and — on Xi and w and M{w) on 2 2 . In the other two problems the "data" of w are 

defined analogously. 

Throughout this paper we assume the boundary 2 of R to be sufficiently smooth so that all 
of the applications of the divergence theorem used are valid. It will be apparent that if the deriva- 
tives of the curvature of 2 are continuous then there is no question of the validity of the indicated 
operations provided the functions themselves are smooth enough. It will also be clear that less 
smoothness of 2 can actually be tolerated in each problem, but the determination of the minimum 
smoothness requirements is not investigated in this paper. 

We make use of the summation convention throughout and employ a comma to denote differ- 
entiation, e.g., w,i=-^-j- We denote the points at which 2i and £2 join as Pi and P 2 . Arc length 
will be measured along 2 in the counterclockwise direction. 

Additional notation will be defined as the need arises in the text. We turn now to the question 
of obtaining bounds in problems I, II, and III. These will be discussed separately in sections 
2, 3, and 4. Certain auxiliary inequalities which will be required in treating the three mixed 
problems will be derived in appendix A. 



2. Problem I 

In a recent paper [4, eq (2.1)] the authors derived the following inequality for a function w 
with piecewise continuous second derivatives: 

f ftt^r ^K^bwHSy +K 2 \[{kw) 2 dAY , (2.1) 

where the constants K\ and K 2 were explicit. (As indicated in [4,5] appropriate constants £1 and 
K 2 are obtained by making use of results of Payne and Weinberger [10]). Using (2.1) as a starting 
point we shall derive the desired L 2 inequality in problem I. 
We first define E(w,w) as 

E(w 9 w)= J [o-(Aw) 2 + (l — (t)w, ijtv, tj]dA. (2.2) 
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It is well known that E(u, u) is proportional to the strain energy. Since (Aw) 2 ^ 2w, ytv, y it follows 
that 



I 



(Aw) 2 dA ^ y^ E(w 9 w). (2.3) 



w)' 



Thus inserting (2.3) into (2.1) we obtain 

( w'cIaY ^ Kil& ufdsl +K 2 Ie(w, «* 

R 1 

where 

k 2 = K 2 [2(l + cr)-T- 

We seek now to bound E{w, w) in terms of the data of problem I. In order to simplify our 
computations we shall assume that 

hm — ^ = 0, lim - , v = (2.5) 

Q-*Pt,Pt dn Q-*PuPs dsdn y ^ J 

where ^ is a point on Si. This restriction is by no means necessary, since bounds could be 

derived which involve the values of t— and : - at Pi and P 2 . In applying the a priori type inequali- 
ty dsdn 

ties, however, w will be set equal to u — 4> where u is the solution to problem I and <£ is an approxi- 
mating function. The function 4> will then be represented as a linear combination of /V linearly 
independent functions <J>, with arbitrary real coefficients. The restriction (2.5) prescribes four 
relations among the coefficients. The Rayleigh-Ritz technique together with these four conditions 
then determine the coefficients. Hence, for practical purposes, (2.5) is no restriction at all. Con- 
dition (2.5) permits us to treat the cases in which Si and £2 do not have end points in common 
(e.g., certain multiply connected configurations) and the case in which points P x and P> do occur, in 
a unified manner. 

In order to bound E{w, w) in terms of the data of problem I, we decompose w as follows: 

w = B + <p (2.6) 



vh 



wUere 



Clearly then 



A 2 £ = 0in R 
B = w on X 

afiffonX, 

1 on2 2 * (2.7) 



AV = A 2 w;inr? 

(f = on 2 

?-0 on I, 

dj\ 

M(<p) = M(w)~ M(B) on 2 2 • (2.8) 
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We now consider the identity 



E(w, w) = E(B, B) + 2E(w, <p) - E(<p, <p) 



and note that 



E(w, <p)=f |£ M(w)dS + [ <pA 2 wdA • 

By Schwarz's inequality we have 

| E(w 9 <p) | ^ { [ (^£)dS [ M 2 (w)dS ] 1/2 + { f ^dA j (A 2 w) 2 dA j 

Since <p vanishes on X, it follows from the duality principle of Fichera [8] that 



1/2 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



where K\ is the same constant as the K\ appearing in (2.1). (In the determination of the explicit 
value for K\ in [5] and [10] the quantity was actually determined for (2.12) rather than (2.1).) It 
is well known also that if <p has piecewise continuous second derivatives in R and vanishes on 
X then 



/.H'H{/. <4 'H" 



(2.13) 



where \\ is the first eigenvalue in the fixed membrane problem for R. If this number is not known, 
lower bounds for it can be easily computed from monotony principles, the Faber-Krahn [7, 9] 
inequality, or other considerations. Inserting (2.12) and (2.13) into (2.11) we obtain 



I E(w, <p) 



[M(w)]*dSy+±{j 



(A 2 w) 2 dA 



1/21 



{J* (AtpfdA 



1/2 



kAJ [M(w)] 2 dsY' 2 +Y {[ (A 2 w) 2 dA 



1/2 



1-r-o- 



E{<p 9 <p) 



1/2 



(2.14) 



the last inequality resulting from (2.3). From (2.14) and (2.9) it follows then that 
E(w,w)^E(B 9 B)+2* 



,1 + a, 

or, using the arithmetic geometric-mean inequality, 

2 



Kx { j [M(w)fdsY 2 ^Y\S {A2w)2dA \ 12 ] [E< *' (p)]1 ' 2 ~ E{(p ' & 



E(w, w) ^ E(B, B) + 



l + o- 



R 



[M(w)fdS \ 2 + 



MJ* 



(A 2 w) 2 dA 



1/2 



(2.15) 



(2.16) 



By combining (2.4) and (2.16) we now have a bound for w 2 dA in terms of L? integrals of the 

J R 

data and E(B, B), where B is the solution to (2.7). We seek now to bound E(B, B) in terms of the data 
of w. 

From the divergence theorem we have 

E(B, B)=\ (ABfdA + (1 - o) \ \B 9 U B, - (AB) 2 ]dA 

J R J R 



-h 



{ABfdA + (l-o) 



dB d 2 B dB d 2 B 1 

ds dsdn dn ds 2 p 



dBV + (§B\ 2 

dn) \ds 



dS. 



(2.17) 
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AT . r dB . dB . v , . , d 2 B ( d 2 B\ . 

IMow it — — and — — are continuous on Z and either - - or — —I is square integrable on z, then we 
ds dn dsdn\ ds 2 J 

may integrate the first (or second) term in the boundary integral on the right by parts on X to obtain, 
for instance, 

We have also made use of the boundary conditions (2.7). An application of the arithmetic-geo- 
metric mean inequality yields for any positive a 

mfl )^£(Ai^ + (l-of£a(^ (2.19) 

An inequality bounding the first term on the right of (2.19) in terms of the data of w is derived in 
appendix A. The insertion of (A. 25) into (2.19) together with (2.7) results in the inequality 

E(B, B) * I £ *B + d, £ (£)"* + <* 3 £ (f* )'« + 1 £ (£)'* + 4 £ (£*)*. (2.20) 

with explicit constants d\. If the region R is simply connected we can actually obtain the inequal- 
ity (see appendix A, in particular (A.26)-(A.28)) 

1S(B,J)<A£ **''***' <& (2.21) 

We now insert (2.20) and (2.16) back into (2.4) to obtain the desired inequality, which may be 
written in the form 

j^dA^A^dS+A^ (g ) 2 dS-M 3 £ (§)'* 



+ 



a[ (I^YdS+^sf (^JjffdS+AJ [M(w;)]^5-fv4 7 [ (\ 2 w) 2 dA. (2.22) 



Here again the constants A\ may be explicitly computed. 

In many of our intermediate inequalities sharper results could have been obtained by leaving 
certain positive weight functions under the integral signs rather than factoring them out. 

It is now possible to employ the technique used in [2] and [3] to obtain pointwise bounds in R 
for w and its derivatives in terms of Li integrals of the data of problem I. It appears likely that by 
using the results of Bramble and Hubbard [1] it would be possible, if sufficient smoothness of the 
boundary and data is assumed, to obtain an inequality of the form (2.22) in which Ai, As, and A& 
are zero. However, the remaining constants would then be more complicated functions of the 
geometry of 2i and £ 2 . 



3. Problem II 

We assume in this section that R and the portion S2 of 2 are such that we can define a vector 
field/' which has piecewise continuous first derivatives in R and satisfies 



/,!-/'/' >P>0 in « 

-fm> g>0on2 2 . (3.1) 
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This is of course always possible if 2 is sufficiently smooth and 2i is nonempty. In fact, if we take 
any sufficiently smooth function <p which is superharmonic in R and whose normal derivative is 
positive on S2 then one such/' is given by 

J <p + a 

where a is chosen so large that <p + a >0 in R. 

In this section we derive a bound for the L 2 integral of an arbitrary C 3 function w in terms of 
Lz integrals of the data of problem II. Again in the application of our inequalities we shall set 
w = u — <$>, where u is the solution to problem II and <t> is the approximating function. In this 
section we assume that O has been so chosen that 

hm UAQ)= hm a v = hm K f • ' = hm v = lim ^ =0. (3.2) 

q^p u p 2 q-+Pi,p> os Q-*p u Pz os Q-^PuPz on Q-*p u Pi oson 



If we choose 

N 

d> = ^a&i (3.3) 

where the 4>, are linearly independent functions, then (3.2) imposes ten conditions on the a,-. If 
2i and 22 have no points in common then conditions (3.2) are unnecessary. 

We now consider the following application of the divergence theorem. 

<bf i n i w 2 dS==\fjw 2 dA + 2 j fww, dA. (3.4) 

Upon rearranging and applying the arithmetic-geometric mean inequality, we obtain 

J \f£-ff)u?dA ^ 6f i mtv 2 dS+ \\&adw\ 2 dA. (3.5) 

Using (3.1) we have 

p J w 2 dA + q I w 2 dS ^ I f i n i w 2 dS J r j |grad w\ 2 dA. (3.6) 

Let us now rewrite (3.4) with w 2 replaced by |grad w\ 2 . It follows then in the same way that 

p |grad w\ 2 dA + q\ |grad w\ 2 dS ^ //i,|grad w\ 2 dS + I w, tjW, X) dA. (3.7) 

J« Ji 2 J^i J R 

Combining the two inequalities and dropping the integrals over 2 2 , we obtain the inequality 

j w 2 dA ^p~ l \ f i n i w 2 dS-\-p- 2 \ /*"/i/|grad w\ 2 dS+ p~ 2 \ w, u w, u dA. (3.8) 

Clearly for cr ^ 

w, ijw, ijdA ^ (1 — o)- l E{w, w). (3.9) 

JR 
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If (r< 0, the factor (1 — a) i must be replaced by 2(2 — a) l . We assume, however, in the following 
that o- ^ 0. Then 

J wHA ^p" 1 ] fmw 2 dS+p- 2 J /V,|grad w\ 2 dS+\p 2 (\ -o-)]~ l E(w, w). (3.10) 



In order to bound E{w, w) in terms of the desired quantities, we decompose w as follows: 



where B satisfies 



w = B -\-(p 

A 2 B = 0'mR 

\w on Si 



(3.11) 



B = < 



on X2 



The function <p then satisfies 



on 2 2 



Afy = A 2 u> in /? 

(£ = — = on Zi 



(3.12) 



(3.13) 



M(<p) = M(w)-M(B), V(<p) = V(w)-V(B)im 2 2 . 

We again make use of the identity 

E(w 9 w)=E{B, B) + 2E(<p, w)-E(<p, if) 



where now 

E(<p,w)=\ j^M(w)dS-( <pV(w)dS+l <ptfwdA. 
We apply the Schwarz inequality as follows: 

E((p, w) 



- I <pAH 

JR 



(3.14) 



(3.15) 



qj i (ffidS+p*f^dA+pqj^ ^dsj 12 ^-^ [M(w)fdS 

+ p- 2 [(tfwfdA+p-'q-i[ (V(w)) 2 ds\ i/2 ' (3.16) 
Now, in view of (3.6), the first term on the right of (3.16) may be bounded by 

q ji (SO dS+ P 2 \/ ldA+ P ( i\ s <P 2 ds^p]Jgrad<p\ 2 dA + qj^ |grad <p\ 2 dS. (3.17) 



dip\ 2 



We have also made use of the fact that It— I ^ |grad </?| 2 . But from (3.7) it follows then that 
pl^ |grad <p\-dA + J |grad <p\ 2 dS *£ I ^if, u <p, ijdA=(l-<T)- l E(<p, <p). (3.18) 



Thus 



|£(<p, w)| ^ {(1 -o-)- 1 ^, ^W^- 1 1 (A/(w)) 2 dS+p- 2 J (tfwfdA+p-iqr 1 ! (V(w)) 2 ds\'\ (3.19) 
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The insertion of (3.19) into (3.14), followed by an application of the arithmetic-geometric mean 
inequality, yields 

E(w, w)^E(B, BHU-o-Hg^f {M(w)) 2 dS+p- 2 f {&w?dA+p- l q- l [ (V(w)) 2 ds\. (3.20) 

A combination of (3.10) and (3.20) thus gives us a bound for I w 2 dA in terms of the data of w and 

JR 

E{B,B). We seek now a bound for E{B,B) in terms of the data of w. 
From (2.17), it follows as before that 



dwA 2 

dn 



dS. (3.21) 



Using the arithmetic-geometric mean inequality we obtain for positive a 8 

e(b, m+fww+a-vy*.^ (igjds+f^-p-*) {^Jds-j^-^Jds. ,3.22) 

As in the previous section we require a bound for I (AB) 2 dA in terms of the data of w. The 

JR 

desired bound was computed in appendix A and is given by (A. 25), i.e., 

{WW * djvdS + d^ (f )'dS + d£ ffidS + ty (HJdS + d^ (f^JdS. (3.23) 

This inequality may be rewritten, using (3.12), as 

fjABWA * Zf*dS + a/ X| (f )'dS + ^ ffidS + ifJ^'dS + &/J*) *■ (3.24) 

A combination of (3.10), (3.20), (3.22), and (3.24) thus leads to the desired inequality, which 
is of the form 

[wHA *s bJ w*dS + ft f (^fJdS + B, f (Jf pS 

W*. i^J dS+B i l {j£> 2dS+B "L 2 lM{w)fdS 

+ B 7 f [V(w)] 2 dS + B s f [tfwfdA, (3.25) 

where the Z?,'s are explicit constants. 

In many cases (3.25) may be simplified by exploiting some of the geometric properties of X. 
For instance, whenever p is positive on Si some quantities may be dropped from (3.21) or used to 
cancel other quantities. 

Pointwise bounds in this pro!' m follow from the results of [2] and [3]. 

4. Problem III 

This problem is the most difficult of the three since for certain geometrical configurations 
it does not have a unique solution. Suppose, for instance, that the domain R lies in the right 
half plane and Si is a portion of the y-axis. In this case the plate may rotate as a rigid body 
about the y-axis. The solution will be unique only up to a rigid body rotation, a term of the form 
w = ax. Note, however, that if 2i is not a straight line then this rigid body motion is no longer 
possible. Because of the inherent difficulty in treating problem 3 we shall derive in this section 
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only an Li bound for the energy integral E(w, w) in terms of the data of the problem. It is possible 
with such a bound, using the techniques of [2, 3], to obtain pointwise bounds for any second deriva- 
tive of w in terms of L 2 integrals of the data. In order to simplify the problem somewhat, we assume 
that w vanishes with its first two tangential derivatives at Pi and Pi. This means merely (as 
pointed out before) that we put certain restrictions on the approximating functions, and it amounts 
to six conditions on the arbitrary parameters in the application of the Rayleigh-Ritz technique. 
With this restriction we decompose w as follows: 



where h satisfies 



w = h + w (4.1) 



A/i = 0in# 

, f w on Si 

h== n v (4-2) 



Then 



& 2 w = A 2 w in R 

w = 0, M(tl) = M(w) - (1 - a) 



2 w , 1 dhl v 



M{w) = M(w) - — - ^ V(w) = V(w) + (l- cr) TTT on S 2 • 
p dn ds 2 dn 

Making use of the triangle inequality, we obtain from (4.1) 

{E(w 9 iv)} 1 * ^ {E(h, h)} 1 ' 2 + {E(w, w)}" 2 . (4.4) 

We compute first a bound for E{h, h). 
Since h is harmonic it follows that 

E(h, h) = (1 - cr) I A, u h, u dA • (4.5) 

J ft 

But, as in (2.17), we observe that 

Applying the arithmetic-geometric mean inequality we obtain 

L A ' <*> »" * J x a 9 (§) 2 dS + j i «*-- p-t) (M )'*-]>' (g)\fi. ,4.7) 

If the region is convex, we could choose a 9 = max p and drop the last two terms. If not, we may use 
the inequality (A. 4) given in appendix A, i.e., 

I (^j 2 dS ^ b i h 2 dS + 6, i (yfdS = K t w 2 dS + b x [ (y J dS. (4.8) 

Upon inserting (4.8) and (4.7) into (4.5), and using the boundary conditions on h, we arrive at the 
inequality 

E{K A) * (1 - <?){ 5o T "^ + 5^ (^) 2 dS + 6 2 f (^Ej ds\ (4.9) 

where 6 , b\ and 62 are explicit. 
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In order to bound E(ib, iv), we decompose w as follows: 

w = B + ip (4.10) 

where 

A 2 B = A 2 w'mR 

M(B)=M(w)onX 

V(B) = l V{iv)on ^ ~ 

' l-F(&) + Ci + C«* + C,yon2 (4.11) 

with Ci, C2, and C ;i chosen to insure existence of the function 5. We indicate now how the C/ are 
to be determined. r r r 

Let us choose the coordinate axes in such a way that I xds = I yds = xyds = 0. Then, 

J 2, J s, J Ii 

in order that (4.11) may have a solution, it is necessary that the following three equilibrium condi- 
tions be satisfied: 

I V(B)dS=[ A 2 BdA 

<f> xV(B)dS-i n. r M(B)dS=\ xA 2 BdA 
Ji Ji Jr 

6 yV{B)dS-i n tJ M(B)dS=l yA 2 BdA. (4.12) 

Ji Js Jr 

These three conditions are now used to determine d, C2, and GY We note first that since h is 
harmonic 



I. 



V(h)dS = 

1 



J xV(h)dS-j n x M(h)dS = Q 
J yF(/i)^S - J n y M(h)dS = 0. (4.13) 

Thus conditions (4.12) are equivalent to 

C,Si = I A 2 zW,4 - f F(w)dS 

C2/1 = J *A 2 wd4 + J n x M(w)dS - [ xV(w)dS 

C3/2 = I yA 2 wGL4 + I n y M(w)dS - J yV(w)dS (4.14) 

where Si denotes the length of Si and 

/i=f ^S, / 2 =f r 2 dS. (4.15) 

Equations (4.14) thus determine the constants Ci, C 2 , and C 3 in terms of the data of problem III. 

In determining the bound for E(w, w) of (4.4), we shall, according to (4.10), need a bound for 
E{B, B). We first indicate how E(B, B) is to be bounded in terms of the data of the problem. We 
present here a method which is applicable only if the domain R is simply connected. The function 
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B is of course determined only up to a rigid body motion, which we assume to be fixe 


d by the 


conditions 

<J> BdS = <t> ^dS = j> y-dS = 0. 


(4.16) 


From the results of [3, eq (4.18)], if R is simply connected, it follows that 




( 2lT ]%ir>r\l ( dB YdS- 2 e<d m 


(4.17) 


{Lj% BdS + % \dn) dS - Dm - a) E{B > B) 


where p 2 is defined in [3, eq (4.31)]. 




From the identity - C C 

6 x i n i B 2 dS = 2\ B 2 dA + 2\ x { BB^dA 




(4.18) 


we obtain, using the arithmetic-geometric mean inequality 




\ B 2 dA ^ r M i B 2 dS + r* \ B, ,fi, x dA • 

J R J 1 J R 


(4.19) 


A repetition of the Inequality with B 2 replaced by j5, ,B, , yields 




f B 2 dA ^ r M <j> B 2 dS + r£ j> | grad B \ 2 dS + r* J tf , y£, u dA • 


(4.20) 


Here nn denotes the maximum value of r= vx i x i on 2. We observed in the proof of (4.17) that 


i»«»*mm'" s 




i | grad B 1 2 dS ^ 4- f B, yB, y<£4 • 

Jl PZJR 


(4.21) 


Thus from (4.20) and (4.21) we have 




1 B 2 dA ^v\B, U B, ijdA ^ -^- E(B, B) 

Jr Jr 1 O" 


(4.22) 


where 




^ = {^(4)'+'i/}(p.)- 1 + r J ' / - 


(4.23) 


By Green's identity 




E(B, B)=\ ^-M{w)dS~ [ BV(w)dS~ [ BV{B)dS+ [ BA 2 wdA= \ ^-M(w)dS 




- f BV(w)dS- I (di + C* + C*y)BdS+ 1 BtfwdA- f ^-M(h)dS+ f flP(A)dS 




- J |^M(w)dS- [ BV(w)dS- f (C, + C 2 * + C 3 y)BdS + [ BA 2 wdA-E(B, h). 


(4.24) 


It follows then from Schwarz's inequality that 




w-t/^-^VH-t/M-^nM-r 




+ J J fi 2 <£4 J ( A 2 w) 2 <//f V 2 + { £(#, B) X* f E(A, A) V'". 


(4.25) 
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Here V(w) is defined as 

\V{w) on 2-2 



V{W) Id + Czx + C 3 y on Si. (4.26) 

From the definition of Ci, C2, and C3 we observe that F(m;) is data. Using (4.17) and (4.22) we 
obtain then 

<«"■ B)) '"' fea^Fi [/«<»>W+(s)T I ,F( "" ) ' <e ])"' 



1-0- J/? 



1/2 



1/2 



(A 2 w) 2 dA + {£(A, A)} ' , (4.27) 



where the bound for the last term is given by (4.9). With (4.27) we now compute the bound for 
E(w, w). 

From (4.10) it follows that 

E(w, w) = E(B, E) + 2E(w, <p) - E(<p, <p) (4.28) 

where 

E(w, if) = I zr M(<p)dS - \ wV(ip)dS + f wA 2 <pdA 
Ji on Ji J/e 

= (4.29) 

from (4.10) and (4.11). Thus 

E(w t w) ^ E(B, B). (4.30) 

Combining (4.4), (4.9), and (4.27) we now obtain (after an application of the arithmetic-geometric 
mean inequality) an inequality of the form 

E(w 9 w)^dJ u?dS + D 2 [ (^jdS + D,[ (^fdS 

+D 4 b [M{w)fdS + dJ> (V(w)) 2 dS + D 6 ( {tfwYdA, (4.31) 



which is the desired result. 



5. Appendix A 



In treating the problems discussed in sections 2, 3, and 4 we required bounds for certain 
integrals of derivatives of biharmonic functions in terms of their Dirichlet data. We derive the 
necessary inequalities here, the most important of which is (A.25). 

Consider the following Dirichlet problem for R: 

A 2 B = 0inR 

B, — given on 2. (A.l) 

on 

Let us compute a bound for {ABfdA in terms of the data of B. We assume that the data are such 

that d 2 B/ds 2 and d 2 B/dsdn are square integrable on 2. 
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Let h be the harmonic function which is equal to B on the boundary. Then from Green's 
identity we have 



f (AB) 2 dA = f [A(B-h)fdA = I ^-(B-h)ABdS. 
Jr Jk Jion 



(A.2) 



Using the arithmetic-geometric mean inequality we obtain, for any positive a , the inequality 



[(hBfdA^i 



(ABfdA < ^ (j> (ABfdS + -?- 
Ir 2. Jx ^«oJi 



dn 



(B-h) 



dS 



^i (ABrds+— 

2 Jx a 



vf) ! - 



dn) 



(A.3) 



In [10, 11] Payne and Weinberger have obtained the following inequality for the last integral 
i (A.3): 

dS (A.4) 



.^rfdS^boi h 2 dS + bi 
i \an/ Js 



a 



jj) 2 ds = bj> B*ds + 



m 



The constants bo and 6i are explicitly determined. We seek now a bound for the first integral on 
the right of (A.3). 

For an arbitrary sufficiently smooth vector field g* and an arbitrary harmonic function H we 
consider the following identity: 



= f g*(B -H), k&BdA = <f «*(fl -H), * £ ( AB)dS 

Jk Jx on 



- g, *(B - H\ fc AS, ufa - g*(fl - //), W A5, t dA. (A.5) 



Integrating by parts further we obtain 

d d 



j^AB 



dx' ' dx k 



(B - H), k dS + l<j> g?n k (AB?dS + f g*(fi - H), * ^~ <% 
■^Jx Jx drc 



+ J g, in k (B - H), iABdS = f { Ag'(B - JSO, i + 2g, K# -H),ik~\g, }Afll Afi<i4 . ( A.6) 

Suppose now that a vector function g* with bounded second derivatives has been found such 
that 

g i = rn (A. 7) 

on 2. II 2eC 2 one possible vector field g l is. given in appendix B. With g l so chosen we take H 
to satisfy 

AH = 0'mR 



^=f-s-4fds 

dn dn Ji on 



(A.8) 



where S denotes the length of 2.. From (A. 7) and (A.8) it follows that the third boundary integral 
on the left of (A.6) vanishes (since B is biharmonic). 
Since g l — m on 2, the first term in (A.6) becomes 



niABl nk 



dx k ) 



ni— l .](B-H), k dS = <pAB 



\ dx 

while the fourth term on the left may be written as 



ds 2 p dn 



(B-H)dS 



m, k n k (B-H)uABdS = i f^ (B-H),iABdS. 

x Jx on 



(A.9) 



(A. 10) 
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Inserting (A.9) and (A. 10) into (A.6), and making use of (A.7) and (A.8), we obtain 



H(AB) 2 dS = <j) AB^- 2 (B-H)dS + S~4 (Afl)p- 1 dS(f ^~dS- 



ds 2 



2 dn 



I on Ids os 



+ nS~4 |5 ds\*BdS + jbB^iB-mti + fyMB- (A.ll) 

We assume now that g l and its second derivatives are uniformly bounded in/?. Then by mak- 
ing use of the arithmetic-geometric mean inequality we obtain for arbitrary positive constants on 



(1 - ai - a 2 - a 3 - ou)i (ABfdS ^ — 



(*3 J2 



fs^-V 



dS + 



(B-H) 

C 2 £ (dB\ 



dS 



-fc(ffds 



ola Ji \dn/ 



Ji \dn/ 
dS + (a 5 + «e + C 3 ) [ ( A£) 2 <£4 



(A. 12) 



+ ^ ((B-H)AB-H),idA+^ \(B-H),ij{B-H),ijdA. 

«5 J/J «6 J/? 

In (A. 12) we have used the abbreviations 

Ci^max -f- — 

C 2 =max -f- m 
z Ion J 

/? 

C 4 =max[A#'A#'] 2 
/? 

C 5 = 4 max|^,ig, |.|. 
/? 

We need now bounds for (B — H), i{B — H), t dA, {B — H),ij{B — H),ijdA and L 2 integrals 

J R J R 

over 2 of the first and second tangential derivatives of H, in terms of the data of B. We consider 

first I (B — H),ij{B — H),ijdA. An integration by parts gives (see 2.17) 
JR 



(A.13) 



j R {B-H) 9ij (B-H),udA = 2i 



d 2 (B - H) d(B - H) 



dsdn 



ds 



dS 



1 



+ j R (ABYdA -<f>^-[(B-H), i(B - H), <]dS. 



(A. 14) 



The first term on the right vanishes in view of (A.8). Thus 
j R (B-H), iA B-H), *U-\WW S-j^dsg^dsJ-^ I , dg 

« f (ABfdA + Cei (jff dS + C 6 



(B-H) 



dS 



TsV-V 



dS (A. 15) 



where Ce is given by 



[Oif p^Oon X 

H f-i) 

P/max 



if p changes sign on 2 
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(A. 16) 



/. 



In a recent paper [6] the authors obtained the following inequality: 

" d 



(B-H),i(B-H\idA^C 7 



dn 



{B-m 



I 



dS+Cs {ABfdA 



(jjffdS + cJ (ABfdA (A.17) 



with explicit constants C 7 and C%. The second inequality results from insertion of the boundary 
condition in (A. 8) followed by an application of Schwarz's inequality. The insertion of (A. 15) and 
(A.17) into (A. 12) then gives 



(l-ai-a 2 -a3-a4)(J> (ABfdS 



1 



aiJs 



ds 2 



(P 2 )max | Ci CgCg , C4C7 



OLi 



a* ae as 



(B-H) 



dB 



-i Ji 



(hi 



dS 



l dS + 



Lr\_ , C5C6 



«3 oc 6 



*<*-« 



dS 



+ 



as a- y 



J. 



We now choose the constants a, in such a way that 

y = I — cx\ — a> — QJ3 — a.| > 0. 



(ABfdA. (A.18) 



(A.19) 



For simplicity we write (A.18) as 



7, d> (A£) 2 dS ^ y 2 



d 2 



WeS"* 



+ y< 



9s 



(fl-tf) 






dS + 7 5 (ABfdA. (A.20) 



We now insert (A.4) and (A.20) back into (A. 3) to obtain 



7i 



2 /J* 



(A£) 2 <£4 ^ y 2 



£<*-«> 



dS+ y.-,+ 



«o/ Jv \3/i/ 



+ 74 



9a 



(B-H) 



*dS + ^$B*dS + ^4 (ffdS. (A.21) 



The constant a () may be taken as 71/7.5 which yields the inequality: 



£■ J (ABfdA ^ 72 <j> 



as 2 



(B-H) 
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2 ^ + (r: i + y 5 )£(|f) 2 ^ 



£<*-*> 



dS + 6075 (J) £ 2 dS + 6,75 



as 



dS. (A.22) 



We need now to bound the tangential derivatives of H in terms of its normal derivatives. 
From the results of [10, 11] we have 

i {f) 2ds * c ° i (f ) 2rfs + ^ I *• * *" * c » £ (If) 2 ^ s ( A - 23) 

with explicit constants C 9 , C 9 , and C 9 . An inequality of type (A.17) has been used to obtain (A. 23). 
We require finally a bound for d) —r dS- 
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Since H is harmonic we obtain, on applying an inequality of Bramble and Hubbard [1], 



d 2 H 

ds 2 



dS ^ Ci 



d 2 H 

dsdn 



dS + Cii 



'dHV 

dn 



dS ^ Ci, 



d 2 B 

dsdn 



dS + C 



i §M 



dBV 



dS (A.24) 



with explicit constants Cio and Cn (see Appendix C). A use of the arithmetic-geometric mean 
inequality in the first and third terms on the right of (A. 22), together with (A.23) and (A.24), yields 
the inequality 



jiABfdA * 2y, j>l^fdS + 2y 2 C l0 j^fdS + (2 y4 + b, y 5 ) £(f J dS 



+ (ys + y,-t2y 2 Cu + 2y 4 C 9 , 



.)£( 



dBV 

dn 



dS+boys&BzdS 



<\>(^JdS + d 4 



dBV 



(A 


.25) 


~dB,idB,~ 


dS 


ds ds 



which was to be derived. 

It might be noted that the left hand side is unaffected by considering B instead of B, where 

B= B + a + aix + a 2 y (A.26) 

In (A.26) the constants a\ are determined by the conditions that 



i BdS= <j> B, x dS= & B, y dS = 



Then, as indicated in [3], it is possible, if R is simply connected, to replace (A. 25) by 



J R 



2 ds ds 



2 ds ds 



(A.27) 



(A.28) 



6. Appendix B. The Vector Field g 

We make the assumption that 2eC 3 and consider the family of parallel surfaces to 2 in R 
given by N(x) = constant, with 2 defined by N(x) — 0. The parallel surfaces need only be defined 
in the neighborhood of £ and we shall be interested in the subdomain of R given by ^ N(x) ^ K~ 1 
where K is a bound on the curvature of 2. The outward normal vector n\ on the parallel surface 
is given by 



m=- 



It is clear then that if we define g l as 

^{(l-K^nuO^N^K- 1 



g 1 - 



then g = ri\ on 2 and g l eC 2 in R. 



otherwise 



7. Appendix C 



The inequality (A.24) follows from eq (2.6) of [1] by making use of [1, eq (2.2)] to prescribe the 
arbitrary vectors g l and symmetric tensors A ij . We sketch here briefly how the argument goes. 

For sufficiently smooth g l and A ij we have the following identity for any harmonic function H: 

j [2{g*A*\ iHtjm-ig-AV), m H, }l ]H, u dA=j> g*AW, iHH^mm-H^m^dS. (C.l) 
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Suppose g l and A ij to be so chosen that on J 

g l = *H 

^ = 6y~^^ (C.2) 

If £ is sufficiently smooth, we may choose g 1 as in appendix B and 

Av^du-gigJ. 

We now decompose the derivatives on the boundary into normal and tangential derivatives to 
obtain 

[n m {8ij — mrij} //, u(2H, jm ni — H, jin m )]dS 

= i \(d 2 H\ 2 (d 2 H\ 2 2dHd 2 H } 1 (dH\ 2 
Jj \\dsdn) \ds 2 ) p ds dsdn p 2 \ds) 

\ ds ds dn dn/ \ ds ds dn dnj J 
Inserting (C.3) back into (C.l), we have upon rearranging and dropping the last two terms 

where £is easily determined once g> and A ij are given. From (A. 8) it follows that 
But from (4.16) 



/. 



'•""-'IS&'-ifflSf+ffi 



dS. 



If we now use the arithmetic-geometric mean inequality in the second term on the right of (C.5) and 
the second term on the right of (C.6), we obtain an inequality of the form 

We now use (A. 23) and apply Schwarz's inequality to the last term on the right after insertion of the 
boundary condition. This leads to (A. 24). 
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